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In this short review presented at FPCP04, Daegu, Korea, we discuss methods leading to 
determinations of a and 7 with practically no theoretical error. The remaining theoretical 
errors due to isospin breaking, neglecting of electroweak penguins or coming from other 
sources are addressed. 



1 Introduction 

We are entering a period of time, when direct determinations of the angles a and 7 of the 
standard unitarity triangle are becoming possible. In this talk we will review the methods 
that are used at present and the related theoretical uncertainties. Surprisingly enough, some 
of the most useful methods were not even talked about before 2003. The questions that will 
be addressed are therefore (i) what is the ultimate precision of different methods and (ii) what 
are we learning about a, 7 now? The last question has been covered in great detail in talks by 
experimental colleagues PP, so only the final results will be given here. 

How can one measure a and 7? The sensitivity to the phases comes from interference. Useful 
methods thus rely on channels with at least two interfering amplitudes and/or interference 
between mixing and decay. In order to extract the weak phases, however, one needs to evaluate 
unknown hadronic parameters that also enter the obsevables. A conservative approach to this 
problem is to extract all the hadronic parameters from experiment. This is accomplished by 
using symmetries of QCD (e.g. C, P, isospin), and by finding channels, where all parameters 
are obtainable from experiment. Another approach is to calculate the hadronic parameters 
using theoretical frameworks like QCD factorization, PQCD, and SCET. This later avenue will 
not be exploited here and the reader is referred to [2] for further details. 

2 Measuring a 

2.1 B(t) -> 7T7T 

This method is due to Gronau and London and dates back almost 15 years ago Let us 
review the method step by step to see where the approximations enter. A completely general 
isospin decomposition of the decay amplitudes is 
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Figure 1: Left: presentation of Eqs. (0), (0), due to Gronau and London [3]. Only one of four 
possible triangle orientations is shown. Right: constraints on a — a e g from isospin analysis [5]. 

where the notation for the reduced matrix elements is Aai- Equivalent relations hold for 5°, 
B~ decay amplitudes A + _, A 00 , A +0 . Note that the AI = 5/2 operators are not present in 
the effective weak Lagrangian, so that A 5 / 2 can only arise from isospin breaking final state 
rescattering effects, such as AI = 2 electromagnetic rescattering of two pions. One can thus 
estimate A 5 / 2 ~ aA\j 2 . Setting A 5 / 2 = therefore means neglecting a ~ 1% correction. Making 
this approximation one obtains two triangle relations 

A + _ + V2A 00 = \ / 2 A +0 , + V2A 00 = V2 A +0 . (2) 

Aside from possible electroweak penguin operator (EWP) contributions, A +Q is a pure tree. 
Neglecting EWP one has an additional relation 

e^A +0 = e-^A +0 \A +0 \ = \A +0 \. (3) 

This allows to extract sin 2a from T(B°(t) — > 7r + 7r~) oc [1 + (7^ cos Ami — SV^ sin Ami] using 
the construction of Gronau and London 0. The observable sin(2a e ff) = S n7T / a/1 — is 
directly related to a through 2a = 2a c $ — 29, where 9 is defined on Fig. 1. 

The difficulty of this approach is the need to distinguish between B — > and B ' -> vr 7r° 
decays, i.e. the need to measure the sides A 00 and A 00 of the triangle relations (see Fig. 1). Since 
the summer of 2004 at least a preliminary isospin analysis is possible, as first measurements 
of B°(B°) -> 7r°7r° rates became available, C o = -0.28 ± 0.39 and Br(B -> 7r°7r°) = (1.51 ± 
0.28) • 10~ 6 [4 . Taking a simple weighted average of Belle and BaBar results on B(t) — > 7r + 7r _ 
(see Table 1), the isospin analysis would at present lead to the constraint on a — a e s shown on 
Fig. 1 [Hj. One does see two emerging peaks when information on Aqq, Aqq is included, however 
more data is needed to constrain a — a e g. At present 

|a-a eff |<39° (90% CL). (4) 

Furthermore the interpretation of a is far from clear due to marginal consistency of S n7T mea- 
surements, Table 1. Recently it was also noted that a e s > a, if the magnitude and phase of 
penguin contributions is not too large j^j. 
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B — > 7T + 7T 


sin 2a c ff 




BABAR [8J 
BELLE 

average 


-0.30 ±0.17 
-1.00 ±0.22 
-0.61 ±0.14 


-0.09 ±0.15 
-0.58 ±0.17 
-0.37 ±0.11 



Table 1: Experimental values of observables in B — > n + n . 

Let us now return to the question of theoretical uncertainties in the isospin analysis [Jj. 
There are two sources of isospin breaking: (i) d and u charges are different, and (ii) m u does 
not equal m<j. The difference of the light quark charges results in additional operators, the elec- 
troweak penguin operators, in the effective weak hamiltonian. Including EWP in the analysis 
will not affect the separate triangle relations (J2J), but only the additional relation Q, since A 0+ 
no longer receives only tree contributions. Remarkably enough, there exist a relation E] 

rrA7=3/2 3 C 9 ± C W 

^eff, EWP - "2 Cl + C 2 

which makes the inclusion of EWP fairly straightforward. Instead of Eq. (J3J) one has 

e*A+ = e-^A+o, (6) 

where 5 ~ 1.5° [llj. The only assumption that entered this estimate is the dominance of EWP 
operators Qg t w, while no estimate of matrix elements is needed. Note that still \A +0 \ = \A +0 \. 
Deviations from this relations would therefore not test the presence of EWP but only the size 
of the Wilson coefficient suppressed EWP operators Qt S . Note as well, that the same relation 
e ny _ e -i(i+s)ji hoidg a i so f or AI = 3/2 (tree) amplitudes in the pp and pre systems. 

Nonzero m u — difference results in tt° — rj — rj' mixing, i.e. 7r° wave function has small 77,77' 
admixtures. Because of this Gronau-London triangle relations (J2J) no longer hold ^2]- Gardner 
[T2"] found that this typically leads to Aa ~ 5° shift in the extracted value of a. Since we now 
have more experimental data about the im system it would be interesting to reevaluate this 
effect, especially if the analysis is extended beyond factorization that was used in ^2]- The 
analysis of 12J also showed that Aa depends on the value of P/T and will thus be different for 
pp, pn systems, where no such quantitative analysis exists at present. 

2.2 Measurement of a from B — ► pp 

The isospin analysis in B — > pp follows the same lines as for B — ► tttt, but with three separate 
isospin relations (J2J), one for each polarization. However, longitudinally polarized final state 
dominates the other two, /+~ = 0.99 ± 0.03+°;°^ [H| and f£° = 0.97±°;°? ± 0.04 HU. This 
simplifies the analysis as there is effectively only one isospin relation. Another difference from 
the 7T7T system is that p resonances have a nonnegligible decay width. The invariant mass 
measured from the decay products can thus differ from the pole mass of the p resonance. The 
two p resonances in the final state can therefore also form an / = 1 state, if the respective 
invariant masses are different ^H]- This affects the analysis at 0(r 2 p /m 2 p ). As shown in ^3] it is 
possible to constrain this effect experimentally by making different fits to the mass distributions. 
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An ingredient that makes the pp system favorable over inr is a small penguin pollution, as 
can be inferred from the bound Br(B -> p°p°) < 1.1 • 1(T 6 (90% CL) 16J (cf. also Fig. 1). This 
gives a measurement of a from S p + P - using isospin analysis [IE] 

a= [96±10±4±11]°, (7) 

with the last error representing the ambiguity due to the presence of penguins. In obtaining 
the above result isospin breaking effects, EWP, non-resonant and 1=1 contributions were 
neglected. 

2.3 B -» P 7T 

Since p^ir^ are not CP eigenstates, extracting a from this system is more complicated. There 
are essentially two approaches, (i) either one uses the full B(t) — ► 7r + 7r~7r Dalitz plot together 
with isospin [T7], or (ii) one uses only the p ± n T region together with SU(3) related modes [T8] . 
If the full B — > 7T + 7r~7r° Dalitz plot is used, one needs to model the Dalitz plot, for instance 
fit to a sum of Breit-Wkrner forms 

f(B° -> 3vr) = BW{s+) A{B° -> p+TT) +W(s_) A(5° -> p-vr + ) +W(s ) -> pV), 

S v ' V «, ' " v ' 

A+ A_ A 

(8) 

where for simplicity only p resonances were kept in the sum, but other resonance can be added. 
From time dependent B(t) — ► 37r Dalitz plot analysis one has 27 real observables, of which 
18 measure the interference between different p resonance bands. In this way it is possible to 
determine A^q, A± t o, up to an overall phase, i.e. there are 11 independent measurables. A 
potential problem can arise from the fact that the peaks of p resonance bands do not overlap, 
but are separated by approximately one decay width. To measure the 11 observables correctly 
one therefore has to model the tails of the resonances correctly. 

In order to extract a from A± ) q, A± <0 additional input is needed. First let us define tree 
and penguin contributions according to whether or not they contain CKM weak phase 

At = e*t± + p±, Ao = e*t + p , (9) 

and similarly for A±,Aq, but with a sign of 7 flipped. The AI = 3/2 part of the weak 
hamiltonian has a CKM phase, so the penguins p±$ are purely AI = 1/2 (neglecting EWP). 
This leads to an isospin relation f7| EI] 

Po = -^{p++P-), (10) 

which reduces the number of unknowns to 10. One possible choice of unknowns is a, \t±\, \to\ 
argt±, \p±\, argp±. There is thus enough information to determine all of them. Explicitly, the 
observable that gives a directly is 

A + _ + A_ + + 2Ax> = T => - Im (|f ) = sin(2«). (11) 

There are some further comments that apply to the Snyder-Quinn method. As already stated, 
the effects due to isospin breaking have not been analysed quantitatively yet. However, isospin 
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breaking will enter only in relation ([10)1 . Since penguins are small, \p±/t±\ ~ 20% ^HJ , it is 
reasonable to expect that isospin breaking effects will also be small, or at least smaller than 
in the B — ► titt case. In addition, if A5/2 7^ 0, only the part of A5/2 that has the same weak 
phase asp± } o will affect the analysis by modifying the relation in Eq. (fTUj) . These contributions 
would come from electromagnetic final state rescatering of penguin contributions, leading to 
negligible effect. 

BaBar performed the Snyder-Quinn analysis (but keeping 10 out of 27 observables fixed to 
zero), obtaining [20] 

a=(113j$±6) Q . (12) 

Note that there is only one solution in [0°, 180°). 

The potential problem of having to model the tails of the p bands can be avoided by using 
just the p ± Tr T final state and the 577(3) related modes [TS] . As in ©, the tree and penguin 
contributions are defined according to their weak phases. In total there are 8 unknowns: \t±\, 
\p±\, arg (p±/t±) , arg(t_/£+), a, but just 6 observables. Additional infromation on penguin 
contributions can be obtained from SU(3) related AS = 1 modes, in which penguins are CKM 
enhanced and tree terms CKM suppressed compared to the p ± 7i T final state. Since penguin 
contributions are small, the error introduced because of the SU(3) breaking will not be large. 
Note that in order to relate the AS" = 1 and AS = channels, annihilation like topologies were 
neglected. 

To resolve ambiguties an additional assumption of arg(t_/t + ) being smaller than 90° had 
to be used. This leads to 

a = (94±4±15)° (13) 

with the last error the combined error coming from a — a e s difference and the estimate of SU(3) 
breaking effects. To obtain this number no interference information was used (i.e. experimental 
data from both BaBar [2II| and Belle |2I] was used). Also, only bounds on penguins were used, 
not a complete SU(3) fit. In the future an unconstrained fit to obtain a could be performed. 
This would lead to a single solution for a, with all ambiguities resolved. As already stated, the 
SU(3) breaking on extracted a would be small, of order p±/t±. A Monte Carlo study with up 
to 30% SU(3) breaking on penguins for instance gives y/{(a ont - a in ) 2 ) ~ 2° [18 . 

3 Measuring 7 

3.1 B ± DK* 

There are many methods that fall into this class, all of which use the interference between 
b — > cus and b — > ucs [22] • In the case of charged B decays this means that the interference is 
between B~ — ► DK~ followed by D — > / decay and B~ — > DK~ followed by D — > /, where / 
is any common final state of D and D. What makes this method very powerful is that there 
are no penguin contributions and therefore almost no theoretical uncertainties, with all the 
hadronic unknowns in principle obtainable from experiment (with problems in measuring color 
suppressed B~ —>■ D°K~ decay [23]). 

Different methods can be grouped according to the choice of the final state /, which can be 

(i) a 

CP- eigenstate (e.g. J\~ 5 7r ) [22J, (ii) a flavor state (e.g. K + ir ) (23) , (hi) a singly Cabibbo 
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suppressed (e.g. K* + K~) 24 or (iv) a many-body final state (e.g. Ks^n^) [25J. There are 
also other extensions: many body B final states (e.g. B + — > DK + ir°) [2S], D°* in addition to 
D°, self tagging D ** [22] or neutral B decays (time dependent and time- integrated) can be 
used (SHI EH!- 

In this talk we focus on extracting 7 from B ± — > (Ks , k + 'k~)dK ± , since this is experimentally 
most advanced. Both experiments use D* and D decays, where a subtelty of a sign flip in the 
use of D* has been pointed out only recently [3U]- The BaBar result [31] 

7 = (88 ±41 ±19 ±10)°, (14) 

should therefore be treated as preliminary only. Belle on the other hand obtains [32"] 

7 = (68+^ ±13 ±11)°. (15) 

Note that only a single solution for 7 is obtained in [0, 180)° range. 

For details on how the method works see [23 EH] E21 EH]- We will just make several state- 
ments regarding the remaining theoretical errors. First of all, it is possible to extend this ap- 
proach beyond Breit-Wigner fits of Dalitz plot, so that there is no modeling error left [23] 133]. 
Also, the effect of D — D mixing is included automaticaly, if D* tagged D decays are used to 
measure the observables of D system. 1 The largest remaining theoretical error is due to possible 
direct CP violation in the D decay, which is, however, highly CKM suppressed by A 5 ~ 5 • 10~ 4 . 
The measurement of 7 will therefore be dominated by experimental errors for years to come. 



3.2 7 from B{t) -»- L>(*)+L>«- 

This is a very recent method [33] . Again, the amplitude is split into tree and penguin according 
to CKM 

A D = A(B° -> D+ir) = ^ + pe^ (16) 

vsy* v; b v vd 

Value of t is determined from B° — ► D^D~ using SU(3) with leading SU(3) breaking correction 
accounted for % = rr-^-, with subleading corrections estimated to be below 5% — 10%. At 

t Vcd JD ' 

present additional approximations are needed to obtain bounds on 7 from D* + D~, D* + D*~ . 
This leads to three viable regions for 7 at 68% CL, 7 e [19.4°, 80.6°] or 7 e [120°, 147°] or 
7 G [160°, 174°]. 



3.3 sin(2/3 + 7 ) 

The combination sin (2/5 ± 7) can (at least in principle) be extracted very cleanly from the 
time dependent measurement B(t) — > D^n~ [3H]. Until the small direct CP asymmetry is 
measured, however, the weak phase 7 and the strong phase 6 can be extracted from S D w±^ = 
2r/(l ± r 2 ) ■ sin(2/5 ± 7 ± 5) only, if the ratio of the two interfering amplitudes r = |v4(I? — > 

D^ +, k~) / A{B — > £)(*) + 7r~)| is known. This ratio can be obtained using SU(3) from Br(B° — > 
1 I thank T. Gershon for pointing this out. 
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Dg + 7t ). Assuming factorization, taking /„(*) / f D (*) from lattice, and neglecting (very) small 
anihhilation like diagrams, this gives = 0.019 ±0.04, rp.^ = 0.017^q;qq7. Using this number 
BaBar obtains | sin(2/3 + t)| > 0.58 (90% CL) from partially reconstructed B -> D^ir* [H7|. 

4 Conclusions 

In conclusion, we have working tools to determine angles a and 7 of the CKM unitarity triangle. 
The experimental situation looks much more favorable than expected a few years ago. For 
instance, measurements of a are already reaching precision level, where one has to start worrying 
about theoretical errors. 
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